AD-A069  991  MASSACHUSETTS  INST  OF  TECH  CAMBRIDGE  LAB  FOR  INFORMA”— ETC  F/G  12/2 
DECENTRALIZED  MAXIMUM  FLOW  ALGORITHMS# (U) 

MAY  79  A SE6ALL  N00014-75-C-U83 

UNCLASSIFIED  LIDS-P-915  NL 


1 OF  1 

AD 

AO 699 91 

a 

Q — - 

jfir  ^ 

IflB 

5WA069991 


nurum  UUV_ 


.1.  REPORT  NUMBER 


«.  TITLE  (and  Sublltla) 


ir.UtUllllNI 


BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.I  S.  RECIPIENT'S  CATALOG  NUMBER 


DECENTRALIZED  MAXIMUM  FLOW  ALGORITHMS 


•».  AUTHOR^ 


Adrian  Segall 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Massachusetts  Institute  of  Technology 
Laboratory  for  Information  and  Decision  Systems  V 
Cambridge,  MA  02139 


tl.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Defense  Advanced  Research  Projects  Agency  . 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 


S.  TYPE  O ? PEPORT  A period  COVERCO 
Paper 

I 

s.  performing  org.  report  number 
LIDS-P-915  ^ 

S.  CONTRACT  OR  grant  NUMBSRf*)  . 

^Contract  ONR/N00014-75-C-lli  3 


13.  PROGRAM  ELEMENT.  PROJECT.  TASK 

AREA  A WORK  UNIT  NUMBERS 

Program  Code. No.  5T10 
ONR  Identifying  No. 049-383 


12.  REPORT  CATE 

May  1979 


IS.  HUMBER  of  pages 
31  pages 


17.  DISTRIBUTION  STATEMENT  (at  tba  mbotro e!  an  farad  In  Block  20.  It  dlltarant  kox  Report) 


ia.  supplementary  notes 


ImSShnH 


Submitted  for  publication  to  "Networks". 


19.  KEY  WOROS  (Continue  on  elde  If  n»e*n«y  end  Identify  by  block  nusibeg) 


20.  ABSTRACT  (Continue  on  geeeree  elde  It  neceeeery  end  Identity  by  btock  nusrbes) 

^^.This  paper  presents  three  algorithms  for  obtaining  maximum  flow  in  a network  using 
distributed  computation.  Each  node  in  the  network  has  memory  .and  processing  capa- 
bilities and  coordinates  the  algorithm  with  its  neighbors  using  control  messages. 

Each  of  the  last  two  versions  requires  more  sophistication  at  the  nodes  than  the 
previous  one,  but  on  the  other  hand  employs  less  communication. 

I* 


DD  IJAN7J  1473/^  EOITIONOF  I NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wtan  Onto  E n farad) 


4 rjL.  fe 





DECENTRALIZED  ^MAXIMUM  JTLOW  ALGORITHMS 


NdmW5-C-1183t  N(WT0l4-77-C-jl532 


Research  supported  in  part  by  the  Defense  Advanced  Projects  Agency 
(monitored  by  ONR)  under  contract  ONR/N00014-75-C-1183  and  in  part 
by  the  Office  of  Naval  Research  under  contract  ONR/N00014-77-C-0532 


V 7 

a 

/ -J— 

t 

1 


» 


f 


i 


1 4 

A 

■J 


Submitted  for 
publication  to 
"Networks" 


DECENTRALIZED  MAXIMUM  FLOW  ALGORITHMS 


Adrian  Seg&ll 

Faculty  of  Electrical  Engineering 
Technion  - Israel  Institute  of  Technology 
Haifa,  Israel. 


Abstract 


The  paper  presents  three  algorithms  for  obtaining  maximum 
flov  in  a network  using  distributed  computation.  Each  node  in  the 


network  has  memory  and  processing  capabilities  and  coordinates  the 


algorithm  with  its  neighbors  using  control  messages.  Each  of  the 


last  two  versions  requires  more  sophistication  at  the  nodes  than 


the  previous  one , but  on  the  other  hand  employs  less  communication 


This  work  was  performed  on  a consulting  agreement  with  the  Laboratory  for 
Information  and  Decision  Systems  at  MIT,  Cambridge,  Mass.,  and  was  supported 
in  part  by  the  Advanced  Research  Project  Agency  of  the  US  Department  of 
Defense  (monitored  by  ONR)  under  Contract  No.  N00014-75-C-1133,  and  in  part 
by  the  Office  of  Naval  Research  under  Contract  No.  0NR/N0001U-77-C-0532. 


1.  INTRODUCTION 


In  this  paper  we  treat  the  problem  of  decentralized  regulation  of 
the  flow  of  a commodity  through  a capacitated  network.  We  assume  that  a 
controller  with  a certain  computation  capability  is  located  at  each  node  of 
a given  network;  each  controller  can  measure  the  amount  of  flow  incoming  to 
and  outgoing  from  the  corresponding  node,  on  each  of  the  adjacent  links  and 
neighboring  controllers  are  able  to  exchange  information  over  the  link  con- 
necting them  in  the  form  of  control  messages.  The  problem  is  to  design  a 
protocol  for  each  of  the  controllers,  so  that  the  combined  algorithm  will 
maximize  the  total  flow  through  the  network  , from  a given  node,  called 
source,  to  another  given  node,  called  sink. 

The  problem  of  maximizing  network  flow  using  central  computation  has 
been  widely  studied  in  the  literature  and  many  algorithms  have  been  designed 
for  its  solution  [1]-  [5].  On  the  other  hand,  technological  developments  of 
mini  and  micro  computers  have  made  it  possible  to  introduce  relatively 
sophisticated  computation  and  large  memory  capabilities  in  each  of  the  nodes 
of  a transportation  or  communication  network.  Decentralized  computation 
provides  a serious  advantage  over  the  centralized  one,  mainly  because  of 
enhanced  survivability  and  reliability,  and  because  it  saves  transmission 
requirements  of  status  information  and  of  command  information  towards  and 
from  the  central  controller.  As  such,  developing  efficient  and  reliable 
decentralized  algorithms  to  perform  various  network  functions  is  of  major 
importance  for  network  design.  The  research  on  decentralized  network  algor- 
ithms is  Just  at  its  starting  stages,  but  several  such  algorithms  have  already 
been  developed  and  validated.  Among  these  we  may  mention  here  algorithms 
for  "HrHmmn  delay  routing  without  and  with  changing  topology  [6]  - [9], 
shortest  path  [10]  - [11]  and  minimum  spanning  tree  [12]  - (1M . 
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Consider  a network  (N,E)  where  N is  the  set  of  nodes,  ECN*N 
is  the  set  of  directed  edges  and  assume  that  the  network  has  no  self-loops 
and  no  parallel  edges.  Two  nodes  s,  teN  are  specified,  where  s is  the 
source  and  t the  sink.  Without  los3  of  generality  we  may  assume  that 
(j,s)  cA  and  (t,j)  i A for  all  J eN.  To  each  edge  (i,k)  eE  ve  attach 
an  integer-valued  capacity  c^.  If  the  capacities  are  not  integers  but 
commensurable,  an  appropriate  scale  factor  will  provide  an  equivalent  problem 
with  integer-valued  capacities.  If  (i,j)  eE  or  (j,i)eE,  then  node  i 
is  said  to  be  adjacent  to  node  j and  the  corresponding  edge  is  3aid  to  be 
Incident  to  node  i. 


A feasible  flow  f is  an  assignment  of  an  integer  f^.  to  each  edge 
( i , j ) such  that 


0 <_  f <_  c for  all  (i,J)  e E and  such  that  for  all 

X J 1 J 


IeN-{s,t}  we  have 


V - l fu  - 0 

J:(i,j)eA  J : ( J ,i )eA 

Each  node  i has  certain  memory,  processing  and  transmission  capabilities 
and  is  assumed  to  know  and  to  be  able  to  update  the  flows  f and  f on 

1 J J i 

all  edges  incident  to'it.  In  addition,  if  (i,j)eE,  then  nodes  i and  J 
are  assumed  to  be  able  to  send  to  each  other  control  messages  that  are  cor- 
rectly received  at  the  other  node  within  arbitrary,  finite,  non-zero  time. 

We  may  note  that  this  assumption  does  not  preclude  that  individual  control 
messages  are  received  In  error  because  of  noise  say.  provided  that  there  i3 
an  error  control  and  retransmission  protocol  that  will  assure  error  detection 
and  appropriate  retransmission  until  the  message  is  correctly  received. 
Observe  also  that  our  assumption  implies  that  if  (i,j)  eE,  then  control 
messages  can  be  sent  over  it  in  both  directions,  even  if  (j,i)  ^ E. 


The  total  flow  F of  a feasible  flow  f is  defined  as  the  net  flow 


from  source  to  sink,  namely 


F - l f * l f 
J : ( 3 , J )cE  9J  J : ( J , t ) cE  Jt 


Cur  goal  is  to  design  a decentralized  protocol  whereby  each  node  i will 

make  computations  and  decisions  based  upon  local  knowledge  of  f and  f 

^ J 

and  upon  control  messages  received  from  adjacent  nodes,  to  achieve  a flow 
for  which  the  total  flow  is  maximum. 


Ji 

f 


In  the  following  sections  we  present  three  decentralized  algorithms 
for  achieving  maximum  flow.  They  are  decentralized  versions  of  the  well-known 
maximum  flow  algorithms  of  Ford-Fulkerson  (FF),  Edmonds-Karp  ( EK } and  Dinic, 
respectively.  All  algorithms  are  based  on  flow  updating  cycles  started  by 
the  sink-node  t,  propagating  in  the  network  while  seeking  a flow  augmenting 
path  from  source  to  sink  and  terminating  back  at  the  sink.  At  the  time  of 
the  termination  of  a cycle,  node  t knows  that  all  nodes  in  the  network  have 
completed  their  part  of  the  previous  cycle  and  are  ready  to  perform  a new 
cycle,  30  that  node  t can  indeed  start  a new  updating  cycle.  The  protocol 
also  Insures  that  at  the  completion  of  the  cycle  that  finally  achieves  maximum 
flow,  the  sink  node  will  be  informed  that  no  more  cycles  are  necessary.  The 
three  algorithms  differ  in  the  procedure  for  finding  augmenting  paths.  In 
the  first  algorithm  (based  on  the  FF  algorithm),  paths  are  searched  and  found 
in  a random  fashion,  employing  a decentralized  protocol  similar  to  the  one 
used  in  [8],  (91  for  adaptive  routing.  In  the  second  algorithm  (based  on  2K) , 
the  augmenting  path  is  always  one  with  smallest  number  of  arcs  and  to  achieve 
this  in  a decentralized  way,  we  employ  simplified  versions  of  algorithms 
proposed  in  [10]  and  [8],  Finally,  in  the  third  algorithm  the  procedure  con- 
sists of  two  kinds  of  actions  as  in  Oinic  [3]:  first  the  nodes  find  the  set 
of  all  possible  paths  with  smallest  number  of  arcs,  then  one  augments  all 


- 5 - 


pc3sible  paths  out  of  this  set,  and  the  procedure  is  repeated  for  longer  and 
longer  paths. 

We  shall  need  now  several  notations  and  definitions.  For  the  purpose 
of  simplicity,  if  (i,J)cE  and  (j,i)^E,  then  we  shall  include  (j,i) 

in  E and  take  Cj^  =0,  fj^  = 0.  For  a given  flow  f,  let  us  define  for 

each  edge  the  quantity 

piJ  = CiJ  " fiJ  + fji  • 

It  can  be  considered  as  the  available  extra-capacity  in  the  direction  from 
i to  J of  the  pair  (i,j),(j,i)  in  the  presence  of  flow  f.  In  particular, 
additional  flow  can  be  pushed  from  i to  J if  and  only  if  p >0.  A 

pair  (i,j),(j,i)  such  that  p > 0 or  p > 0 (or  both)  will  be  called 

J J 1 

a link,  and  if  p > 0,  we  shall  say  that  the  link  is  outgoing  from  i and 
incoming  to  J.  Observe  that  these  definitions  are  flow  dependent  and  that 

t 

if  P^j  > 0 and  pj^  > 0,  then  the  link  is  both  outgoing  from  and  incoming 

to  i (as  well  as  from  and  to  J).  For  a given  node  i,  we  shall  use  the 

notation  [i,j]  for  the  pair  (i,j),(j,i)  and  also: 

ri  ■ <k  i v * 01  ■ 

°i-  °>  • 

Given  any  two  nodes  i,k  in  the  network  N,  an  augmenting  path  from 

i to  k is  defined  as  a series  of  distinct  nodes,  i = i ,i, ,i_,...,i  = k 

o 1 2 m 

such  that  for  all  n * 0,1, . . . ,(m-l)  there  is  a link  outgoing  from  i and 

n 

incoming  to  i , . . The  quantity  p = min  p.  . , where  the  minimum  is 

n+l  i ,1  , 

n n+l 

over  n * 0,  ...  (m-l),  will  be  called  the  available  capacity  of  the  augment- 
ing path.  If  there  is  at  leasx  one  augmenting  path  from  i to  k,  we  say 

that  i has  access  to  k and  k is  accessible  from  i.  We  denote  by  A 

— — — — — — — * r 

the  set  of  nodes  that  have  access  to  sink  at  the  time  of  starting  the 

r-th  cycle,  and  by  L the  set  of  links  connecting  any  two  nodes  in  A at 

r . r 


*■.  . -f  # T ' v m -rm. 


time  t . 
r 


r , 
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In  each  of  the  following  sections,  we  shall  first  describe  the 
general  distributed  protocols  corresponding  to  each  of  the  versions,  then 
indicate  the  exact  algorithm  to  be  performed  by  each  node  as  a part  of  the 
protocol  and  finally  prove  that  the  protocols  indeed  achieve  maximum  flow. 
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3.  VERSION  1 

The  protocol  described,  in  this  section  is  the  decentralized  version 

of  the  Ford-Fulkerson  labeling  procedure.  As  said  before,  the  algorithm 

consists  of  a series  of  cycles  started  by  and  terminated  at  the  sink  node  t. 

We  assume  that  before  the  first  cycle  starts,  a feasible  initial  flow  exists 

in.the-aetwoEkicevgrk.f:,  ■•=  0 for  all  <•  edges  land  also  .all:nodese-,i  iare-,  in 

lk 

state  51  with  the  same  cycle  number,  n^  = 0 for  all  i.  A cycle  consists 
of  two  phases*  Phase  A is  started  at  sink  and  propagates  into  the  network. 
During  this  phase  augmenting  paths  from  network  nodes  to  sink  are  being  built 
in  such  a way  that  all  paths  belong  to  a single  directed  tree  rooted  at  sink. 
All  nodes  that  have  access  to  the  sink  t will  eventually  enter  the  tree, 
and  therefore  the  source  node  s will  also  join  the  tree,  provided  that  the 
existing  flow  is  not  maximal , Phase  B of  a cycle  propagates  over  the  tree, 
from  its  leaves  cowards  the  sink.  During  this  phase,  the  nodes  along  the 
source-sink  augmenting  path  update  the  flows  (or  equivalently  the  quantities 
p ) and  all  other  nodes  make  no  flow  changes.  However,  they  participate  in 
phase  B in  order  to  complete  their  part  of  the  cycle  and  to  finally  inform 
the  sink  that  the  present  cycle  has  been  terminated.  3y  the  time  the  propaga- 
tion of  phase  B arrives  at  the  sink,  all  nodes  in  the  tree  have  had  completed 
their  part  of  the  current  cycle,  and  therefore  the  sink  knows  that  a new  cycle 
can  be  started. 

To  accomplish  its  part  of  the  protocol,  each  node  i keeps  the  fol- 
lowing set  of  variables: 

n^  - cycle  identification  number  (values  0,  1); 

a^  - source-sink  augmenting  path  identification  (values,  0,  1); 

d^  - available  capacity  of  augmenting  path  from  i to  sink 
(strictly  positive  integer -valued,  except  that  d * ®) ; 


;* 


p - pointer  to  nexo  node  on  augmenting  path  from  i to  sink 
(called  father) ; 

N^(k)  - for  each  adjacent  node  k (values  0,  1). 

The  meaning  qf-r^(k).  will  be  explained  presently.  ;The. propagation  , 
of  the  cycle  is  achieved  using  control  messages  sent  by  nodes  to  their 
adja:ent  nodes  and  having  the  format  MSG(n,d,a),  where,  if  i is  the  send- 
ing node,  then  n = n.  , d = d. , a = a.  . The  variables  n and  n.  serve 

i i i 1 

to  distinguish  a given  cycle  from  the  previous  one.  The  variable  a^  takes  on 
normally  the  value  0;  it  changes  to  1 if  node  i’  finds  out  that  it 

* 

belongs  to  the  current  augmenting  path  from  source  to  sink.  If  a.  =0, 
the  quantity  d^  denotes  the  available  capacity  of  an  augmenting  paoh  from 
i to  3ink.  If  a^  * 1.  then  d^  indicates  the  available  capacity  of  the 
augmenting  path  from  source  to  sink  and  it  is  exactly  the  amount  of  flow 
pushed  from  source  to  sink  during  the  current  cycle. 

We  next  proceed  to  describe  the  propagation  of  a flow  increasing 
cycle  in  the  network  and  the  actions  taken  by  the  nodes  to  participate  in 
such  a cycle  (the  numbers  in  parentheses  refer  to  the  Tables) . Just  before 
the  cycle  is  started,  all  nodes  in  the  network  are  in  state  SI  with  identi- 
cal cycle  numbers  {n^}.  A cycle  is  started  by  the  sink  t when  it  per- 
forms the  transition  T12  (cf.  Fig.  l)  from  state  SI  to  state  S2  (1.51).  At 
this  time  it  sets  the  identification  number  of  the  new  cycle  n as  the 
binary  complement  of  the  number  of  the  previous  cycle  (1.52)  and  sends  a 
control  message  on  each  of  its  incoming  links  (1.53).  The  new  cycle  identi- 
fication number  will  be  carried  by  all  control  messages  of  the  present  cycle. 

An  arbitrary  node  i performs  its  part  of  phase  A when  it  receives  the 
first  message  belonging  to  the  present  cycle  (1.13).  At  this  time  it  updates 
its  own  cycle  number  n^  (1.17),  it  chooses  the  originator  of  this  message  t 
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is  father  (l.l8)  and  sets  the  available  capacity  to  sink  on  the  chosen 

path  d.  to  the  minimum  of  d»  (received  in  the  message  as  d)  and  of 

i J 

ohe  available  capacity  p ^ on  link  (i,£).  Then,  to  propagate  phase  A,  it 
sends  a control  message  over  all  incoming  links,  but  not  to  its  father  p^ 
(1.20).  It  is  easy  to  see  that,  as  proven  in  the  Appendix,  all  nodes  having  an 
augmenting- -path' to  sink  will  indeed. perform  T12.and  the  links  {(i,p^)}. 
will  form  a tree  rooted  at  sink. 

Having  completed  its  part  of  phase  A,  node  i waits  for  phase  B,  and 

this  will  be  performed  when  node  i had  received  a control  message  of  this 

cycle  over  all  of  its  incoming  links.  Since  node  i is  not  interested  in 

messages  coming  on  links  that  are  not  incoming,  such  messages  are  immediately 

"bounced"  by  node  i (1.4,  1.5).  To  recognize  the  situation  that  will 

enable  it  to  perform  phase. B,  node  i stores  the  cycle  number  of  any  received 

message  from  t in  IT  (£)  (1,3)  and  will  indeed  perform  phase  B when 

N (k)  = n for  all  incoming  links  (1.21).  Also,  if  in  the  meantime  it 
2 2 

receives  a message  with  a = 1,  node  i recognizes  that  it  belongs  to  the 

* 

source-sink  augmenting  path  for  this  cycle,  performs  the  appropriate  flow 
increase  (1.6  -1.10)  and  during  transition  T21,  which  is  exactly  phase  3, 
it  augments  the  flow  on  the  link  (i,p_*)  (1 .22  - 1 .21*) . In  any  case  it  sends 

iuring  T21  a control  message  to  p^,  thus  informing  it  of  the  completion  of 
its  part  of  the  cycle.  It  is  easy  to  see  (cf.  Appendix)  that  this  informa- 
tion will  propagate  downtree  and  by  the  time  the  sink  nor'e  t performs  T21 
(I.5U),  all  nodes  on  the  tree  will  have  completed  their  part  of  the  present 
cycle.  If  an  augmenting  path  has  been  found  and  the  sink  still  has  incoming 
links,  then  a new  cycle  will  be  started  (1.52);  otherwise,  maxi  mum  flow 
has  been  achieved.  Finally,  observe  that  the  source  node  s performs  the 
same  operations  as  all  other  nodes,  except  that  it  has  no  incoming  links  and 
that  it  is  always  on  the  source-sink  augmenting  path  (1.27-1.42). 


J 
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The  decentralized  protocol  described  above  indeed  achieves  maximum 
flow.  To  prove  this,  it  is  sufficient  to  3hov  that  for  any  cycle  started  by 
the  sink  node  at  time  say; 

- the  cycle  will  be  completed  in  finite  time; 

- if  the  flow  is  not  maximal  at  time  x , the  cycle  will  increase  . 
the  flow  by  a strictly  positive  integer  amount; 

- if  the  flow  is  maximal  at  time  tn  , the  cycle  will  not  change 
the  flow; 

- if  at  the  end  of  the  cycle  the  flow  is  not  maximal,  a new  cycle  will 
be  started  by  sink  at  that  time. 

The  exact  proofs  of  these  properties  are  given  in  the  Appendix. 
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Edmonds  and  Karp  [2]  have  shewn  that  if  we  insist  that  the  chosen 
augmenting  path  from  source  to  sink  has  minimum  number  of  links,  then  the 
number  of  augmentations  is  0 ( j N | ^ ) , as  opposed  to  the  Ford-Fulkerson  algor- 
ithm where  this  number  can  be  bounded  only  by  the  value  of  the  maximal  flow. 

Ve  shall  new  present  the  decentralized  version  of  the  Edmonds-Karp  algorithm. 

As  in  Version  1,  the  present  protocol  consists  of  consecutive  cycles, 
started  by  and  completed  at  the  sink  node  t.  A cycle  cons.sts  of  two  phases, 
the  first  one  forming  the  tree  of  shortest  paths  to  the  sink  and  the  second 
one  propagating  dovntree  to  finally  inform  the  sink  of  the  completion  of  the 
cycle  while  increasing  the  flow  on  the  augmenting  source-sink  path.  Here 
one  finite-state  machine  for  the  node  needs  three  states  (see  Fig.  2):  the 
state  S_.  of  a node  i is  SO  before  it  hears  of  a new  cycle,  it  is  SI 
whenever  the  node  participates  in  the  first  phase  of  a cycle,  namely  while 

it  looks  for  its  shortest  path  to  sink,  and  - = 52  after  the  node  has 

found  thi3  path  and  ’until  it  completes  the  second  phase  of  the  cycle.  This 
version  differs  from  Version  1 mainly  in  the  procedure  for  the  node  to 

choose  its  father  p , namely  during  the  time  the  node  is  in  state  SI.  In 

Version  1,  node  i chooses  as  the  link  over  which  it  receives  the 

firat  message  of  the  present  cycle.  If  the  delays  on  all  links  were  identi- 
cal, then  clearly  Version  1 would  also  provide  the  link  corresponding  to  the 
shortest  path.  Since  thi3  is  not  the  case  however,  the  idea  of  the  procedure 
hare  is  to  have  the  protocol  act  as  if  the  delays  were  identical.  Explicitly, 
during  the  phase  of  looking  for  the  3horcest  path  (i.e.  in  state  SI),  a node 
i sequentially  learns  (via  control  messages  from  its  neighbors)  whether  or 
not  the  sink  i3  at  distance  1,  2,  3,  ate.  from  it.  Suppose,  for  axample, 
that  node  i has  already  learned  that  the  3ink  is  not  at  distance  1 or  2. 

If  one  of  the  nodes,  k 3ay,  such  that  [i,k]  is  outgoing  from  i,  reports 


to  i that  k is  at  distance  2 from  sink,  then  node  i knows  that  its 
shortest  path  to  sink  is  via  k and  is  3 links  long.  The  protocol  provides 
this  information  by  requiring  each  node  i to  keep  the  following  variables 


in  addition  to  n^,  d^,  p_^,  which  have  the  same  purpose  as  in  Version  1: 


d^  = » - if  shortest  path  has  not  been  found  yet;  otherwise  same 

meaning  as  in  Version  1; 

z^  - distance  from  node  i covered  by  the  algorithm  until  now 

(values  0,1,. . . |n| ); 

M^(k)  - for  each  adjacent  node  k,  counting  the  number  of  messages 

received  from  k during  the  present  cycle;  it  is  exactly 

the  last  z^  reported  by  k to  i and  received  at  node  i, 
(values  0,1, . . . , | N | ) ; 

D^(k)  - for  each  adjacent  node  k (values  0,®). 

The  meaning  of  D^(k)  will  be  explained  presently.  The  quantity  M^(k)  is 
initialized  to  (-1)  at  the  end  of  the  previous  cycle  (11.36)  and  is  incre- 
mented (II. U)  every  time  node  i receives  a message  from  k.  The  counter 
z.  is  3et  to  0 when  node  i enters  the  present  cycle  (11.22),  incremented 
to  1 after  having  received  at  least  one  message  from  all  outgoing  links,  i.e. 
after  M^k)  > 0 for  all  3uch  links  (II. 24),  (11.27),  incremented  to  2 after 
M^(k)  >_  1 and  so  on  (11.25),  (11.29).  Every  time  z^  is  incremented,  a 
control  message  i3  sent  over  all  incoming  links,  to  inform  the  neighbors  that 
i has  covered  all  nodes  that  are  z^  or  less  links  away.  Node  i declares 
that  it  has  found  the  link  (i.p^  providing  the  shortest  path  to  sink  when- 
ever it  has  received  from  p^  a message  with  d # « (II. 5).  (II. 6),  (11.27), 

and  it  is  ready  to  increment  z.  to  z +1  * M(p, ) + 1 (11.27).  At  this 

1 Pi  1 

time  it  performs  the  transition  to  state  S2  and  from  here  on  acts  as  in 

« 

Version  1. 
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In  order  to  validate  Version  2,  it  is  sufficient  to  show  that  for 
all  nodes  i that  have  access  to  the  sink  at  the  beginning  of  a cycle,  the 
link  (i,p.)  as  chosen  by  the  algorithm  during  this  cycle,  provided  the 
shortest  augmenting  path  to  the  sink.  This  is  because,  except  for  the  phase 
of  chcosing  (i,p^),  the  algorithm  acts  as  in  Version  1 and  that  version 
would  choose  (i,p^)  as  providing  the  shortest  path  if  control  message 
delays  were  identical  on  all  links.  Therefore,  since  validation  of  Version  1 
va3  proved  for  arbitrary  control  message  delays,  the  case  'under  considera- 
tion here  becomes  a special  case. 

Lemma  1.1 

For  each  node  i with  acces  to  sink,  the  link  [ijP^J  at  the  time 
when  i performs  transition  T12,  corresponds  to  the  augmenting  path  with 
smallest  number  of  links  connecting  i to  sink. 

] 

Proof 

We  proceed  by  induction.  With  the  notations  at  the  end  of  Section  2, 
consider  all  nodes  in  A^  that  have  an  augmenting  path  with  m links  or 
les3  to  3ink.  Suppose  that  these  nodes  satisfy  the  statement  of  the  Lemma. 

Then  a node  i whose  distance  to  sink  is  (m+l)  must  have  an  outgoing  link 
to  a node  k with  distance  m and  no  outgoing  link  to  nodes  with  distance 
less  than  m.  The  link  [i,k]  provides  the  shortest  augmenting  path  from 
i to  sink,  and  when  k performs  T12  it  will  send  its  (m+l)-st  message  to 
node  i with  d ^ This  will  make  d^  f «,  and  when  node  i will  also 
have  received  at  least  (m+l)  messages  from  its  other  outgoing  links,  it 
will  perform  T12  with  p * k. 


5,  VERSION  3 

p 

Dinic ' s algorithm  [3]  achieves  maximum  flow  in  0(|n|  ) augmenta- 
tions by  first  finding  the  set  of  all  shortest  augmenting  paths  from  source 
to  3ink  in  the  subnetwork  (A^L^),  and  then  looking  for  augmenting  paths 
only  in  this  subnetwork.  For  a given  flow,  the  ’union  of  all  shortest  source- 
sink  paths  is  called  the  referent  of  the  network  and  after  exhausting  all 
paths  in  one  referent,  a new  referent  with  longer  paths  will  be  searched  for. 
Suppose  that  for  the  given  flew  under  consideration,  the  shortest  source- 
sink  path  has  m links.  Then  the  referent  contains  (m+l)  levels,  where 
t is  in  level  0,  the  source  s is  in  level  m,  and  an  arbitrary  node 
i e is  in  level  J if  its  shortest  path  to  sink  has  J links,  and  if 
there  is  an  incoming  link  [£,i]  to  i such  that  i is  in  level  (j+1) 
of  the  referent.  In  this  case  link  [£,i]  belongs  to  the  referent  (for 
details  see  [3]). 


As  in  the  previous  sections,  the  distributed  version  of  Dinic's 
algorithm  is  composed  of  a series  of  cycles.  Each  cycle  here  consists  of 
fir3t  finding  the  referent  (Part  A in  Table  3),  and  then  of  a series  of  sub- 
cycles during  which  augmenting  paths  in  this  referent  are  searched  for.  During 
the  first  part,  each  node  finds  out  if  itself  and  which  of  its  adjacent  links 
belong  to  the  referent.  A node  i in  the  referent  compiles  the  sets  Rin 
and  R . of  nodes  that  are  one  level  above  it  and  one  level  below  it, 
respectively.  Only  those  nodes  (and  the  appropriate  links)  will  participate  in 
the  second  part  of  the  cycle,  during  which  30urce-sink  augmenting  paths  belong- 
ing to  the  referent  will  be  obtained  in  the  same  way  as  in  Section  3 of  this 
paper.  (Part  B in  Table  3).  We  therefore  have  to  describe  here  only  the 
first  part  of  a cycle,  i.e.  determining  the  referent,  which  is  done  as  fol- 
lows. First  a node  will  find  its  distance  from  sink  in  the  same  way  as  in 
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Section  U,  so  that  in  state  SO  and  SI  it  will  behave  as  in  Table  2.  The  dif- 
ference is  in  state  S2,  because  we  require  that  a node  will  positively  know 
whether  it  belongs  to  the  refer ent  or  not  at  the  tine  of  performing  T20,  and 

if  it  dce3,  which  of  its  incoming  and  outgoing  links  are  in  R and  R , 

in  out 

respectively.  If  a node  i is  at  distance  J from  the  sink,  then  it  will 
belong  no  level  J of  the  referent  if  indeed  it  enters  the  referent.  In 
this  case,  the  adjacent  links  that  will  enter  the  referent  will  be  outgoing 
lizuu  to  the  nodes  distanced  (j-l)  hops  from  sink  and  incoming  links  from 
nodes  distanced  (j+1)  hops  from  sink.  Observe  that  in  Section  U,  if  D^(k)  * 0, 
namely  if  a message  with  d + • has  been  previously  received  from  k,  then 
this  means  that  k has  found  its  distance  from  sink,  which  is  exactly  equal 
to  M^(k).  New,  node  i finds  out  its  distance  from  sink  in  transition  T12. 
Therefore  any  ke  0^  such  that  D^(k)  * 0,  M^k)  = z^  (here  we  mean  the 
value  of  z^  Just  before  the  transition)  will  positively  enter  the  referent 
if  i does  and  all  links  [i,k]  will  also  enter  the  referent  in  the  direc- 
tion from  i to  k.  We  indicate  this  in  (III. 3).  In  addition,  if  such  a node 
k is  also  in  I , node  i will  not  send  a message  to  it  while  performing  T12 
vis  in  Version  2),  but  rather  when  performing  T20  (III. 5),  (III. 6).  3y  that  time, 
i will  know  whether  it  has  entered  the  referent . Node  i will  actually  enter 
the  referent  by  command  received  over  an  incoming  link  (III. 2)  and  performing 
T20  without  receiving  such  a command  is  the  information  that  node  i did  not 
enter  the  referent.  The  variables  a and  a play  here  the  same  role  as  in 
Sections  3,*S  except  that  here  they  indicate  that  node  i belongs  to  the  refer- 
ent, whereas  previously  they  indicated  that  i is  on  the  chosen  source-sink  path. 


From  the  description  of  the  algorithm  it  is  clear  that  Part  A of  a cycle 
of  the  distributed  algorithm  here  will  give  the  same  result  as  the  construction 
l of  ,a  referent  in  [3],  and  this  validates  our  algorithm.  Details  of  the  proof 

are  essentially  identical  to  the  proofs  in  the  Appendix  and  therefore  deleted. 


L 
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6.  Summary 


We  have  described  distributed  versions  of  three  well-known  maximal 
flow  algorithms.  Probably  other  existing  algorithms  [U] , [5]  can  be 
decentralized  using  similar  ideas.  As  mentioned  before,  the  processing 
requirements  are  simpler  in  Version  1,  become  mere  complicated  in  Version  2 
and  even  more  so  in  Version  3.  On  the  other  hand,  as  presently  seen,  the 
communication  requirements  are  reduced. 

For  integer  valued  capacities.  Version  1 needs  a maximum  of  F 

max 

augmentations,  where  F is  the  value  of  the  maximal  flow.  Since  the  algor- 

max 

itba  requires  a maximum  of  1 message/edge  (in  each  direction)  for  each  augmenta- 
tion, the  communication  requirement  is  0(F  ) messages/edge.  Version  2 

max 

requires  0 ( | N | ^ ) augmentations,  and  since  the  maximum  distance  to  sink  is 
jN|,  each  augmentation  requires  no  more  than  |n|  +1  messages /edge.  Hence 
the  communication  requirement  is  no  more  than  0 ( | N | ^ ) messages  edge.  Finally 
in  Dinic's  algorithm  [3]  the  number  of  required  cycles  is  no  more  than  | N [ . 

In  order  to  find  the  referent  one  needs  no  more  than  | JT | +1  messages/edge 
as  in  Version  2,  and  the  number  of  possible  subcycles  is  bounded  by  |sj  , 
since  each  subcycle  saturates  at  least  one  link.  For  each  3ubcycle  one  needs 
no  more  than  1 message/edge  as  in  Version  1,  and  hence  the  communication 
requirements  are  bounded  by  0( |n|  • ( |s|  + | N j ) ) = 0 ( | N | | E | ) <_  0 ( | N | ^ ) 
messages /edge. 
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In  this  Appendix  we  shall  prove  that  the  protocol  proposed  in  Sec- 
tion 3 indeed  achieves  maximum  flow.  The  other  Versions  have  been  validated 
in  the  body  of  the  paper.  We  shall  first  need  several  definitions  ana  nota- 
tions. Recall  that  a cycle  is  started  by  the  sink  t when  it  performs  T12 

and  is  completed  when  3lnk  performs  T21.  Let  x^  be  the  time  when  the  r-th 

cycle  is  started  and  the  time  of  its  completion,  which  is  the  same 

as  the  time  of  starting  the  (r+l)-st  cycle  (see  1.56).  We  say  that  the  net- 
work (N,E)  is  0-relaxed  at  time  x if  all  nodes  icN  are  in  state  SI, 

r 

if  no  messages  are  in  transit  and  if  all  nodes  i in  Ap  have  cycle  number 
n^  * 0.  Similarly  for  1-relaxed.  Also,  in  order  to  specify  the  value  of  a 
variable  at  a given  time  r,  we  put  t in  parentheses,  e.g.  n^(x)  denotes 
the  value  of  n^  at  time  r. 

Suppose  that  the  network  is  0-relaxed  at  time  x^  and  a cycle  is 
started  by  sink  at  this  time,  i.e.  sink  performs  T12  and  changes  a from 

w 

0 tc  1.  In  order  to  see  exactly  what  happens  within  a cycle,  let  us  assume 

that  whenever  (1-5*0  will  hold  (if  at  all),  instead  of  starting  a new  cycle 

after  performing  T21  as  required  by  (1.56),  the  sink  will  stop.  We  shall 

show  that  indeed  (1.5*0  will  hold  in  finite  time  after  t , at  time  x 

r r+i 

say,  and  at  that  time  the  network  will  be  1-relaxed.  We  shall  also  see  the 

operations  performed  by  the  nodes  during  the  period  from  x^  to  x^+^  and 

that , unless  the  flow  was  maximal  at  time  x , the  cycle  will  have  strictly 

r 

increased  the  flow  by  an  integer  amount.  Then,  by  induction  on  the  cycles, 
m*Y-i rmim  flow  will  indeed  be  reached  within  a finite  number  of  cycles. 

Lemma  1 

Suppose  that  at  time  x^  the  network  is  0-relaxed  and  sink  performs 
T12  (while  changing  n from  0 to  1).  Suppose  that  if  after  time  x the  sink 

X It 
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perf-.rms  T21,  it  will  not  start  a r.ev  cycle,  but  will  ranher-  stop.  Then, 
with  the  notations  of  Section  2, 


ii)  each  node  i must  change  n^  from  0 to  1 before  it  generates  any 

message; 

(b)  each  node  icA  will  perform  T12  at  most  once;  nodes  not  in  A 

r r 

perform  no  transition; 

( c ) all  generated  messages  carry  n * 1; 

(d)  each  node  ie  A will  perform  T21  at  most  once; 

r 

(e)  no  more  than  one  message  is  sent  on  each  link  of  L in  each  direc- 

r 

tion.  No  messages  are  sent  on  links  not  in  L . 

r 

Proof 

Observe  that  by  the  assumption  of  the  Lemma,  at  time  t all  nodes 

r 

i£ Af  have  n^  * 0,  sc  that  (a)  is  clear  from  (1.13),  (I. IT).  Now  we  prove 

(b)  - (c)  by  a common  induction.  For.  any  instant  x > x , suppose  that 

r 

(b)  - (c)  hold  until  time  T-.  We  shall  show  that  whatever  happens  at  time  x 
will  not  violate  (b)  -(c).  Suppose  the  first  part  of  (b)  is  violated  at 
time  x,  i.e.  node  i performs  T12  for  the  second  time.  Then  (1.13)  says 
that  it  receives  at  this  time  a message  with  n * 0,  which  must  have  been 
generated  before  time  x-,  violating  (b)  before  x-.  Since  messages  are 
sent  by  i either  to  k e 1^  (1.20),  or  on  links  from  which  messages  have 
been  received  ((1.5)  and  (1.25) ).  nodes  not  in  A^  receive  no  messages  and 
hence  perform  no  transitions,  completing  the  proof  of  (b).  Now  suppose  node 
i generates  a message  with  n ■ 0 at  time  x,  violating  (c).  Then  it  must 
have  n^(x)  * 0 violating  (b)  before  or  at  time  x. 

Nov  (d)  is  clear,  3ince  T12  and  T21  must  alternate  (Fig.  1).  In 
order  to  prove  the  first  part  of  (e),  we  suppose  again  that  it  holds  through- 
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cut  the  networlc  until  time  r-  and  show  that  at  time  t it  cannot  be  violated. 
First  observe  thac  if  condition  (I.U)  holds  at  time  x-,  then  by  (c)  we  have 
n(x-)  * n^(x-)  ■ 1,  so  that  i must  have  performed  T12  earlier.  Therefore 
if  (I.U)  holds,  then  £ # p^,  since  otherwise  p^  must  have  sent  two  messages 
to  i violating  (e)  before  x- . Now  node  i sends  a message  to  p^  in  T21, 
a message  to  each  k e I , k ? p,  in  T12  and  a message  to  each  'mi  I , k f p. 
in  (1.5) > so  that  the  first  part  of  (e)  holds.  The  second  part  of  (e)  follows 
from  the  second  cart  of  (b). 

□ 

Lemma  2 

Suppose  the  same  assumptions  as  in  Lemma  1 hold.  Then 

(a.)  each  of  the  nodes  in  A and  only  those  nodes  will  perform  T12  exactly 

r 

cnce; 

(b)  the  set  and  the  links  (i,p^)  form  a tree  spanning  (A^jL^)  and 

rooted  at  the  3 ink  node  t ; 

(c)  each  node  is  A will  perform  T21  exactly  once; 

r 

(d)  the  last  node  to  perform  T2I  is  the  sink  and  at  that  time,  xr+j_  say, 

the  network  i3  1 -relaxed  and  A , CA  . 

r+1  — r 

Proof 

(a)  At  time  x all  nodes  i in  A are  in  state  SI  with  n,  = 0.  At 
r r i 

time  sink  performs  T12,  changes  n^  from  0 to  1 and  sends  a message 

with  n » 1 over  incoming  links.  Therefore  the  nodes  i at  the  other  end 

of  these  link3  will  eventually  receive  a message  wixh  n * 1,  and  unless  they 

have  performed  T12  before,  they  will  perform  T12  now  (1.13)  while  changing 

n^  from  0 to  1.  Each  such  node  J will  send  a message  with  n * 1 on  all 

of  its  incoming  links,  except  to  p which  ha3  already  performed  T12,  so 

J 

that  by  induction,  all  nodes  that  have  at  time  x^  at  least  one  augmenting 


-~T.T,r 

V . 
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path  to  sink  will  perform  T12  and  change  n^  from  0 to  1.  The  above  reason- 
ing also  shows  that  only  nodes  in  Af  will  perform  T12. 

(b)  The  graph  formed  by  the  links  (i,p.  ) clearly  spans  (A  ,L  ).  If 

i r r 

L is  chosen  by  i as  p^  as  in  (1.18),  then  L must  have  performed  T12 
before  i,  and  hence  thi3  graph  is  loop-free.  Since  each  i has  a unique 
p^»  the  graph  is  a tree.  Since  there  must  be  at  least  one  node  i such 
that  p * t , and  t has  no  "father"  p , the  tree  is  rooted  at  t . 

(c)  Consider  first  a node  i which  is  a leaf  of  the  tree  (i.e.^k  s.t. 
p,  3 i).  Let  us  look  at  an  arbitrary  node  k such  that  ksl  , k 5*  p., 

£ X X 

to  which  i sends  a message  when  performing  T12  (1.20).  When  this  message 
arrives  at  k,  it  must  hold  that  n = n^,  otherwise  k will  choose  i as 
p violating  the  fact  that  i is  a leaf.  Therefore,  if  k^I  then  (1.4) 

£ X / 

will  hold  and  k will  send  a message  to  i.  On  the  other  hand,  if  ke I. , 
then  k has  sent  a message  to  i while  performing  T12.  Therefore  i will 
eventually  receive  messages  from  ail  ks  I , (1.21)  will  hold,  i will 

perform  T21  and  send  a message  to  p . Similarly,  by  induction,  we  can  show 
that  each  node  in  A^  will  perform  T21,  this  action  propagating  down-tree. 

(d)  The  reasoning  in  (c)  shows  that  the  3 ink  t i3  the  last  to  perform 

T21.  If  a node  i^A  , then  it  and  all  nodes  to  which  it  has  access  (which 

r 

are  also  cot  in  Ap)  performed  no  transitions,  hence  remained  in  state  SI 

and  had  no  change  of  flow,  implying  that  they  are  not  in  A either. 

r+i 

Finally,  if  ieA^,  then  it  has  received  all  messages  on  links  [i,k]  such 

that  k e I before  performing  T12.  If  k^I^,  then  i e 1^  and  k has 

received  a message  from  i before  k performed  T21.  But  i can  send  such 

a message  only  when  receiving  a message  from  k.  Therefore  i receives  all 

messages  sent  to  it  before  sink  performs  T21,  so  that  there  are  no  messages 

in  transit  at  time  T and  hence  the  network  is  1-relaxed. 

r+1  n 
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Theorem 

Let  F be  the  total  flow  entering  the  sink.  In  the  protocol  indi- 
cated in  Section  3 and  Table  1,  we  have  for  all  cycles 

F(tr+1)  >F(xr)  (A.l)  i 

with  equality  if  and  only  if  F(t  ) is  maximal.  If  the  latter  holds,  then 

r 

the  sink  will  stop  after  the  (r+l)-st  cycle  or  possibly  even  before  this 
cycle.  Therefore  maximal  flow  will  be  reached  in  a finite  number  of  cycles. 

Proof 

Clearly  the  flew  is  maximal  at  x if  and  only  if  s i A . If 

r r 

seA  , it  will  eventually  enter  the  tree,  terform  .T12  and  immediately  T21 
r 

while  sending  MSG(n  ,d  ,a  =1)  to  o (i.Ul).  When  k = p receives  it, 

s s s s s 

it  sets  d,  = d and  a,  = 1 and  when  performing  T21  sends  MSG(n,  ,d,  ,a,  =1) 

& s k k k k 

to  p . By  induction,  the  nodes  of  the  entire  source-sink  path  on  the  tree 
& 

will  perform  similar  operations,  while  also  changing  the  flow  (or  equivalently 

the  available  extra-capacities  p ^) . Therefore  (A.l)  holds  with  inequality. 

. 

On  the  other  hand  if  s^A  , either  I at  time  x is  empty  in 

r Xi  x* 

which  case  the  sink  starts  no  new  cycle  or  else  the  cycle  is  started,  but  s 
does  not  enter  the  tree  and  the  cycle  will  be  completed  with  a = 0.  In 
eitner  case  the  flow  is  maximal  and  the  sink  stops  triggering  cycles  (I. 57). 

□ 
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£ 

if 


IA3LS  1:  Algorithm  for  Version  1 


la) 

Algorithm  for  Node  isN-{s,t} 

Note:  Message  MSG(n,d,a)  received  from  node  t is 
algorithm  in  the  form  MSG(n,d,a,£) . 

1.1 

Message  Handling 

1.2 

When  receiving  MSG(n,d,a,£) , execute: 

1.3 

N^-d)  * n; 

I.U 

if  n = n . , til . ( comment : t t n . , a 

= o), 

1.5 

then  send  MSG(n. ,d. ,0)  to  d: 

1.6 

if  a = 1 (comment:  i is  in  state  S2. 

t i p ) 
i 

1.7 

then  o .£  *■  o + d. 

1.8 

P-di  °ti  ~ d> 

1.9 

\ - d. 

1.10 

ai  "1; 

1.11 

execute  Finite-State-Machine. 

1.12 

Finite-State-Machine  Transitions  (see  Fig. 

1) 

J 


Hote:  The  Finite-State-Machine  is  executed  until  no  more  transitions 
are  possible. 


1.13  T12 

Condition  12 

MSG(n  # ,d,a,d) . 

I.lU 

Comment  12 

a " 0.  at  = 0,  ou  > 0. 

1.15 

Action  12 

h 4 

I.l6 

°i 4 * 0>; 

1.17 

ni  n; 

1.18 

Pt  ♦ t.\ 

1.19 

di  f min(d,pi^) 

1.20 

send  MSG(n1>d1 .a^)  to 

ke  Iit  k j p^. 
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I * 21 

T21 

Condition  21 

Vkel.,  -chen  N.(k)  = n . 
i ii 

1,22 

Action  21 

If  a.  = 1,  then 
l 

1.23 

1. 24 

1.25 

V*V,dl’ 

°iPi  ~ °iPi  " di’ 
send  MSG(n^ ,d^ ,a^)  to 

1.26 

a.  *■  0. 
l 

lb) 

Algorithm  for  Source  Node  s 

Note 

: Xs  = * 

1.27 

Message  Handling 

1.28 

When 

receiving  MSG(n,d,a,£) , execute: 

1.29 

if  n = n , 

3 

1.30 

then  send  MSG(n  ,d  ,0;  to  t\ 
s s 

1.31 

execute  Finite- 

-State-Machine  for  s. 

1.32 

Finite-State-Machine  Transitions  for  s (see 

1.33 

T12 

Condition  12 

MSG(n  * n^ ,d,a,£) . 

1.34 

Action  12 

d min(d,p  , ) ; 

S SJL 

1.35 

n *■  n; 

S 

1.36 

P3  - u 

1.37 

a *■  1; 

S 

1,38 

execute  transition  T21. 

1.39 

T21 

Condition  21 

none. 

1. 40 

1.41 

Action  21 

P «-  P - d ; 

sp  sp  s 

3 3 

send  MSG(n  ,d  ,a  ) to  p 

q e e i 

a *-  0. 
s 


I.U2 
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TABLE 

2:  Algorithm  for  Version  2 

Algori 

thm  for  Mode  i s N - (s  ,t} 

II.  1 

Message  Handling 

II. 2 

When  receiving  MSG(n,d,a,£) , execute: 

i 

i 

II. 3 

if  d f 

® , then  D.  it)  «-  0, 

II,!* 

M (Z)  - 
i 

M (l) + 1, 

II.  5 

if  d * 

®,  d = ®,  then  p.  +•  t,  d.  min(d,p.»}, 

II.  6 

if  d i 

®,  d^  # ®,  M^(£)  < M(pi),'  then  p^  *■  t,  d^ 

min{d,p  , 

11,7 

if  n = 

n^ , i i I , d p « , i is  in  state  S2  or  SO 

(comment : t ? p . , a 

= 0); 

II. 8 

then 

send  MSG(aJ,d.,o)  to  l\ 

II.  5 

if  a = 

1 (comment:  i is  in  state  S2,  & ? p.) 

II.  10 

then 

V‘“l£‘d' 

II.il 

p4l-»£i-d. 

11.12 

"*•  dL. 

11.13 

ai  " 1; 

II.  iU 

execute 

Finite-State-Machine 

11.15 

Finite-State-Machine  Transitions  (see  Fig.  2) 

11.16 

TOi  Condition  Cl  MSG(n  t n ,d,a,£). 

11,17 

Action  01  I - (k|p^  > 0}; 

11,18 

°i  ^ {k^pik  > °}l 

11.19 

ai  ♦ 0; 

11.20 

n; 

11.21 

d *■  ®; 

11.22 

zi  «■  0; 

11.23 

send  MSG(ni,4^,a^)  to  all  kel^. 

4 





# 

■gap  « — < 

■ 
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II  24 

Til 

Condition 

11 

{Vl-ceO,,  then  M.(k)  > z>  and 

i i — i 

{[d^  = «]  or  [(d^-)  and  ( < M_.  (p^ )] } 

11.25 

Action  11 

* 1; 

11.26 

send  MSG(n^ ,d = «,a^)  to  all  kel  . 

II.  27 

T12 

Condition 

12 

{t^keO,  , then  M,(k)  > z.}  and 

i i i 

{d_.  # »>  and  ^ = M (p  )}. 

11,28 

Comment  12 

a . =0. 
l 

11,29 

Action  12 

V11 

11,30 

send  MSG(n^ ,d^ ,a^)  to  all  ke  1^,  k / p . 

11.31 

T20 

Condition 

20 

l/kel.  , then  D (k)  = 0. 

11,32 

Action  20 

i i 

if  a.  = 1,  then 

— l 

11,33 

II.  3U 

V'V**1' 

11.35 

send  MSG(n.  ,d  ,a  ) to  p_,  ; 

11.36 

fksl  1^0  , then  M (k)  «•  -1,  D(k)  •; 

Note:  To  save  3pace,  we  do  not  indicate  explicitly  the  algorithms 
for  source  and  sink.  The  source  s acts  as  all  other 
nodes,  except  that  it  has  no  incoming  links,  so  that  it 
performs  T20  Immediately  after  T12  and  also  ag  = 1.  The 
sink  t has  only  transitions  T02  and  T20  which  are  identical 
to  T12  and  T21  of  Version  1.  It  will  stop  the  algorithm 
when  it  will  have  a..  * 0 at  the  time  of  performing  T20. 

Then  maximum  flew  is  achieved. 


TA3LE  3:  Algorithm  for  Version  3 


Part  A:  Finding  the  Referent. 

All  instructions  are  exactly  as  in  Table  2,  except  for  the  following 
changes : 


III.l 
III.  2 


In  (II.  7)  change  L 1 I to  £jfIUR  . 

i i out 

Change  (II. 9)  - (11.13)  to:  if.  a = 1 (comment : i is  in  state  S2) 

then  a.  1 and  include  L in  R.  . 

i in 


Action  12  should  read: 

III. 3 RQut  = (k|keO.,  D (k)  = 0,  M.(k)  = * };’ 

III.U  2i  zi  + 

III. 5 send  MSG(n,  ,d.  ,a.  ,r  ) to  all  ke  I,  , k^R  . 

1 1 I 1 I • out 

III.  6 In  Action  20,  delete  (H.  32)  - (II.  3*0  and  change  (11.35)  to: 

send  MSG(n  ,d. ,a. ) to  all  'k  e R . . 

i i l out 

III. 7  In  Action  20,  add:  n^  *■  0. 


Part  B:  Finding  Augmenting  Paths 

The  algorithm  is  almost  identical  to  Table  1.  All  variables,  states 
and  transitions  here  will  have  the  same  names  as  in  Table  1,  except 
that  they  will  be  primed  (e.g.  n^ , si',  T1'2',  etc.)  to  distinguish 
them  from  those  of  Part  A.  All  instructions  are  exactly  as  in 
Table  1,  except  for  the  following  changes: 

111. 8 In  (1. 4)  change  Li  I to  eR^. 

111. 9 Change  (1.8)  to:  «■  P^  - d"  and’ then  if  p^  = 0,  'exclude  L 

from,-  Rifl, 

111. 10  Delete  (1.15),  (I.l6)  from  Action  12. 


III. 11 

Change 

(1.20) 

to: 

send  MSG' (n^ ,d^,a^ ) 

to  all  keR,  . 

in 

III. 12 

Change 

(1.21) 

to: 

VkeR.  , then  N’(k) 
’ in  i 

* ni' 

III. 13 

Change 

^ * 

(1.24) 

to: 

P .■  D.  - d and 

ipt  ipt  i 

then  if  p.  = 0, 

ipi 

l : 

from  R 

• 

out 

III. 14 

Change 

(1.57) 

to: 

else  perform  T01. 

T21 


Fig.  2 - Finite-State-Machine  for  Version  2. 


Til 


exclude 


Fig-  3 


Finite-State-Machine  for  Version  3. 
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